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With high-Tc cuprates in mind, properties of correlated d x 2_ y 2-wecve superconducting (SC) and antifer- 
romagnetic (AF) states are studied for the Hubbard (t-t'-U) model on square lattices, using a variational 
Monte Carlo method. We employ simple trial wave functions including only crucial parameters, such as 
a doublon-holon binding factor indispensable to describe correlated SC and normal states as doped Mott 
insulators. U/t, t' /t and 8 (doping rate) dependence of relevant quantities are systematically calculated. As 
U /t increases, a sharp crossover of SC properties occurs at U co /t ~ 10 from a conventional BCS type to a 
kinetic-energy-driven type for any t' /t. As S decreases, U co /t is smoothly connected to the Mott transition 
point at half filling. For U/t ^ 5, steady superconductivity corresponding to the cuprates is not found, 
whereas the d-wave SC correlation function P^° rapidly increases for U/t ^ 6 and becomes maximum at 
U — U co - Comparing the <5 dependence of P/j° with experimentally observed dome-shaped T c and conden- 
sation energy, we find that the effective value of U for the cuprates should be larger than the band width, 
for which the t-J model is valid. Analyzing the kinetic energy, we reveal that for U > U co only doped holes 
(electrons) become charge carriers, which will make a small Fermi surface (hole pocket), but for U < U co 
all the electrons (holes) contribute to conduction and will make an ordinary large Fermi surface, which is 
contradictory to the feature of cuprates. By introducing a proper negative (positive) t'/t, the SC (AF) state 
is stabilized. In the underdoped regime, the strength of SC for U > U co is determined by two factors, i.e., 
the AF spin correlation, which creates singlet pairs (pseudogap), and the charge mobility dominated by 
Mott physics. In this connection, we argue that the electrons near the antinodal points in the momentum 
space play a leading role to stabilize the d-wave state, contrary to the dichotomy of electronic roles in the 
momentum space proposed for the two-gap problem. We also show the instability of the hole-doped AF 
state against phase separation. 

KEYWORDS: superconductivity, cuprate, antiferromagnetism, pseudogap, doped Mott insulator, doublon-holon 
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1. Introduction 

The essence of superconductivity (SC) in the high-T c 
cuprates is most probably captured by two-dimensional 
(2D) t-J and Hubbard models; 1 ' 2 ' a variety of theo- 
ries have been developed on the basis of these mod- 
els. 3 ~ 5 ' In the ground state of the f- J- type models, 
which deal with strongly correlated regimes, different 
approaches concluded that superconducting (SC) states 
with d x 2_ y 2--wave symmetry are stabilized in the param- 
eter range relevant to the cuprates, as we will see further 
in §2.1. Meanwhile, approaches starting from the weak- 
correlation limit 6, 7 ' of the Hubbard model provided com- 
plementary, but often qualitative, views of SC for the 
cuprates, particularly in thermal and dynamical proper- 
ties. Nevertheless, quantum Monte Carlo (QMC) meth- 
ods, 8 ~ 10 ' quantitatively reliable approaches for finite cor- 
relation strengths, unanimously yielded a result that the 
SC correlation function never develops in weakly corre- 
lated regimes as compared with the noninteracting case, 
and therefore cast a serious doubt on the ground of the 
weak-correlation theories for applying them to realistic 



correlation strengths. Thus, it is important to clarify how 
properties of a d 2 .2_ J; 2-wave SC state evolve as the corre- 
lation strength is increased in the Hubbard model. 

To this end, the variational Monte Carlo (VMC) 
method 11,12 - 1 is useful for its continuous applicability 
from weakly to strongly correlated regimes in a quan- 
titative manner. This is by virtue of correct treatment 
of the local electron correlations. Since the early stage of 
research in the cuprates, VMC studies 13-15 ' have mainly 
treated i-J-type models and yielded results basically con- 
sistent with the behavior of cuprates. 5 ' Later, the present 
authors have studied a many-body SC state, applying 
it to the 2D Hubbard (t-t'-U) model [U: onsite Coul- 
omb correlation, t, t'\ hopping integrals to the nearest- 
neighbor (NN) and diagonal sites]. 16 ~ 19 ' In refs. 16 (for 
t' = 0) and 17, we discussed superconductor (SC)- 
insulator (Mott) transitions at half filling arising at U = 
U c ~ 7t without directly introducing an AF correlation. 
According to these and recent studies, 20,21 ' a binding ef- 
fect between a doubly occupied site (doublon, D) and an 
empty site (holon, H) plays a leading role in nonmag- 
netic Mott transitions. Ref. 16 also showed for doped 
cases that a, sha,rr> crossover exists at U — IL„ ~ W 
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[W(= 8t): band width], where the nature of the d x 2_ y 2- 
wave SC state, ^q, changes qualitatively from a usual 
BCS type for U < U co to an unconventional strongly 
correlated type for U > U co , in which is stabilized 
by the reduction in kinetic energy. This crossover for SC 
is smoothly connected to the Mott transition at half fill- 
ing (U co — > U c ), as the doping rate 5 is decreased. Thus, 
the strongly correlated SC state for U > U co deserves a 
doped Mott insulators, 22 ) and will be deeply influenced 
by the effect of D-H binding correlation. 

In this paper, as a development of refs. 16 and 17, we 
study SC and antiferromagnetism (AF) independently 
for doped cases by systematically changing the values of 
U/t, t'/t and 5, using D-H binding wave functions. The 
primary purpose is to establish that the nature of SC 
is qualitatively different between the regimes of U < U co 
and U > U co by studying various quantities, and that the 
behavior of cuprates is contradictory to that of d-wave 
BCS-type SC realized for U < U co in various aspects, but 
consistent with the behavior of doped Mott insulators 
(U > U co ). 

The second purpose of this paper is to consider a dif- 
ferentiation or dichotomy of electronic roles in the mo- 
mentum space, 23 ' 24 ) in connection with the pseudogap 
problem. Early experimental studies in the so-called two- 
gap problem 25-28 ) posed a question on a familiar view 
that the pseudogap stems from incoherent pairs of elec- 
trons 29 ~ 32 ) for T c < T < T*. Two characteristic en- 
ergy scales (gaps) were discovered whose (5-dependences 
arc mutually opposite in the underdoped regime. They 
seemed to stem from the electrons with wave numbers 
near the nodal (tt/2,tt/2) and antinodal (n, 0) points, 
and correspond to the SC gap and pseudogap, respec- 
tively, indicating that the origin of the pseudogap is not 
directly connected to pairing. Recent experiments for the 
pseudogap have taken on more complicated aspects, 33 ~ 4 °) 
and are still actively developing. 41, 42 > In this study, we 
argue that the two energy scales are interpreted as a 
pairing gap and its correction by the charge fluctua- 
tion released from a Mott insulating state. And that the 
electronic states near the antinodal point (it, 0) [not in 
the nodal region <~ (tt/2,-jt/2)] are crucial to the cZ-wave 

sc 32,43) 

A part of the results were reported previously. 44,45 ) 
The rest of this paper is organized as follows: In §2, 
we explain the model and method used in this paper. 
In §3, we present the results for t' = to grasp basic 
features of U /t and 5 dependence of relevant quantities, 
and argue the existence of a crossover and its relevance to 
cuprates. In §4, we discuss the effects of diagonal transfer 
t', which provides a clue for understanding the stability 
of SC and AF states. In §5, we argue that the electrons 
near the antinodal point play a leading role to realize 
the d-wave SC. In §6, the main results are summarized. 
In Appendix A, we give supplementary explanations of 
the model. In Appendix B, we compare some forms of 
doublon-holon binding factors as a supplement to §2.2. 
In Appendix C, we describe in detail the definition and 
justification of the long-distance value of d-wave pair cor- 
relation function in the present cases. 



2. Formalism 

In §2.1, we explain the Hubbard model and the back- 
ground of this study. In §2.2, we introduce many-body 
trial wave functions used in this paper. In §2.3, we briefly 
summarize the setting of VMC calculations. 

2.1 Hubbard Model for Cuprates 

As a model of the CuC-2 planes, we consider the Hub- 
bard model on a square lattice: 

H = H kin + Hu 

kcr j 

with dj = rij^riji (doublon projector) and rij a = c^ a Cj a . 
Although the band structure e(k) is qualitatively im- 
portant, 46 ~ 52 ) we consider, for clarity, only the minimum 
terms to distinguish material dependence i.e., NN (t) and 
diagonal-neighbor (t 1 ) hoppings: 

e(k) = — 2i(cos k x + cos k y ) — At' cos k x cos k y . (2) 

We use t and the lattice constant as the units of energy 
and length, respectively. As discussed in Appendix A, 
the effective onsite repulsion U/t may differ between 
electron-doped and hole-doped cases. For the ease of cal- 
culations, we can restrict the electron density to n — 
N/N s < 1 (N: electron number; N s : site number), or the 
doping rate to (5=1 — n > by a canonical transfor- 
mation. Then, t'/t becomes negative (positive) for hole- 
(electron-)doped cases (see Appendix A). 

The studies of this model for the cuprates are roughly 
classified into weak- and strong-correlation theories. In 
the former, 6 ^ stability and some properties of d x ?_ y i- 
wave SC are discussed, using RPA, 53 ) FLEX approxima- 
tions, 54 -* renormalization-group 7 ) and perturbative 55,56 ) 
approaches. As we will show in §3 and §4, however, no 
substantial enhancement of the SC correlation function 
appropriate for the high-T c cuprates is found in a weak- 
correlation regime (U/t ^ 5). Such a result is not re- 
stricted to our studies. Early QMC studies 8 -* drew neg- 
ative conclusions for the appearance of SC in the weak- 
coupling regime (U/t = 2-4). A recent QMC study 9 ) us- 
ing an expansion with Gaussian bases also obtained a 
negative result for SC in the overdoped regime up to 
U/t = 7 for t' = 0. Another recent study 10 ) concluded 
that the d x 2_ y 2-wave pair susceptibility is not enhanced 
for t' = 0, but a SC Kosteritz-Thouless transition ex- 
ists for moderately large U/Vs and t'/t = —0.2. Thus, 
it is improbable that the sm&ll-U/t Hubbard model with 
t' = gives such a sufficiently high T c as the cuprates 
display; a doubt arise as to the reliability of the weak- 
correlation theories for the cuprates, at least in a quan- 
titative sense. In this connection, we will check whether 
or not the introduction of t'/t term can enhance SC in 
the weak correlation regime. 

In the strong-correlation limit, the Hubbard model is 
mapped to t- J- type models, in which robust d x 2_ y 2-wave 
SC is found by VMC, 13-15 ) exact diagonalization, 57 ) a 
projector Monte Carlo methods, 58 ) besides by various 
mean-field-type approximations. 4 ' 5 ) In the correspond- 
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ing large-[//i regime of the Hubbard model, firm d- 
wave SC and its competition with AF orders were found 
near half filling by VMC studies 59-61 ) using trial func- 
tions similar to ours. 16-19 -* Another promising quanti- 
tative method to treat this regime is extensions of the 
dynamical mean field theory (DMFT). First, using a dy- 
namical cluster approximation, a finite temperature in- 
stability to d-wave SC was found. 62 ' Subsequently, the 
stability of (I- wave SC and the competition with AF were 
studied using different variations of DMFT. 63 " 65 ) As far 
as we know, there are no theories that reliably drew neg- 
ative results for SC in this regime. It seems appropriate 
to assume that robust d x 2_ y 2-wave SC corresponding to 
the high-T c cuprates is realized in the Hubbard model 
with U ^ W near half filling. 

Finally, we touch on the coexistence of <i-wave SC and 
AF long-range orders. This problem has been repeatedly 
addressed for the 2D t-J and Hubbard models, using 
e.g. VMC methods 59 ' 66 ' 67 ) and DMFT. 63 ' 68 ) These stud- 
ies found the stability of coexisting states against both 
pure d-wave and AF states near half filling; other stud- 
ies 65, 69, 70 ) showed that such mixing states appear for 
small S only under certain conditions. Experimentally, 
explicit evidence of microscopic coexistence has not been 
discovered for most high-T c cuprates, 71 ) except for multi- 
layered Hg systems. 72 ~ 74 ) It is probable that the coexis- 
tence can be realized in clean Cu02 plains in the multi- 
layered systems, while the absence of the coexistence in 
usual cases will be due to disorder. Meanwhile, inhomo- 
geneous phases such as stripe 75 - 1 can compete or coexist 
with the d x 2_ y 2 wave SC. 41,42,76 ) Anyway, as a basis to 
tackle complicated states, it is vital to study the pure 
(i-wave and AF states separately and their competition. 

2.2 Variational wave functions 

We implement a series of VMC calculations. Remov- 
ing irrelevant variational parameters, we elucidate basic 
behaviors of simple correlated d-wave SC and AF wave 
functions of the Jastrow type, 77 ) 



(3) 



when it is applied to the Hubbard model, eq. (1). Here, V 
denotes a many-body correlation (Jastrow) factor com- 
posed of projection operators, and $ a one-body (mean- 
field) wave function given by a Slater determinant. 

In the many-body part, we consider two projection op- 
erators: V — VqVq- 16 ' 17 ^ Although the onsite repulsive 
(Gutzwiller) projector, 78 ) 



V G = l[[l-(l-g)d j } 



(4) 



with < g < 1 is of primary importance, intersite corre- 
lation factors are indispensable for a qualitatively correct 
description of Hubbard-type models. 79 ) Among them, an 
attractive factor between a doublon and a holon is crucial 
at half filling, which is defined as, 80 ' 81 ) 



Vq = ]J(1- Qj) • 



(5) 



where Q 



(=Q?) 



is a D-H projector of the site j: 



diHO—h j+ r) + hj]l(l-d j+T ) 



(6) 



Here, hj = (I — rij^)(l — nj^), and r runs over all the 
NN sites of the site j. The projector Vq with Q| yields 
(1 - ^) jVd+Mi , if isolated Af B doublons and Af H holons 
exist in the electron configuration, where an isolated dou- 
blon (holon) indicates a doublon (holon) not accompa- 
nied by holons (doublons) in its four NN sites. The pa- 
rameter n (< 1) controls the strength of D-H correla- 
tion; for /j = 1, a doublon and a holon are completely 
bound in mutually NN sites, for fi = 0, they are com- 
pletely free, and for /j < 0, they becomes repulsive to 
one another. As we repeatedly showed, Vq or its ana- 
log plays a leading role to induce conductor(metal)-to- 
nonmagnetic-insulator (Mott) transitions 17 ' 20 ' 21 ' 82 ) and 
spin-gap transitions in attractive Hubbard models. 83 ' 84 ) 
Recent studies using more elaborate long-range D-H fac- 
tors 20 ' 21 ) showed that the simplest NN form works 
unexpectedly well. We use this form Q| for 5 = in 
this study. We found that a D-H correlation factor is 
still important for 5 ^ 0, 16 ) if 5 is roughly within the 
range of SC for the cuprates. A NN D-H pair yields, 
with a single hopping process, an antiparallel-spin pair, 
which contributes to the pairing in <i-wave SC. Because 
the symmetry between a doublon and a holon is broken 
for S 7^ 0, an asymmetric form of Qj seems reasonable. 
Thus, for S > 0, we use a simple asymmetric form: 



Qi = Qf =V dj - h j+7 



(7) 



for the ease of calculations. In Appendix B, we check 
different forms of Qj in detail. The main conclusion is 
that the details of Qj makes no difference at least qual- 
itatively; a simple form of Qj like eq. (7) preserves the 
essence of the D-H binding mechanism. In §3.1, we will 
discuss the condition that the D-H binding factor be- 
comes significant. 

For the frustrated cases (f ^ 0), D-H correlation be- 
tween diagonal sites may also play a certain role. There- 
fore, we similarly introduce a D-H projector Qj with a 
variational parameter \x' for inter-diagonal sites. Eventu- 
ally, the D-H factor for t' ^ becomes 



^Q=Il( 1 -^)( 1 -^) 



(8) 



In most calculations for t' — 0, we fixed fi' at zero, be- 
cause the optimized value of /j,' is small and the energy 
gain thereof is very small (~ 10~ 4 i even for large U /i). 85 ) 
Now, we turn to the one-body function $ in eq. (3). As 
a continuation of the preceding papers, we mainly study 
a fixed-A d^-^-wave singlet (BCS) state, $d(Ad, C): 86 ) 



|0>, 



with 



a k = — = 



rr. 



(9) 



(10) 
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where £ is a variational parameter which is reduced to the 
chemical potential for U/t — > 0. Since it is known that 
the d x 2_ y 2-wave is the most stable among various gap 
shapes near half filling, 13 ' 14 ) here we exclusively treat 
the basic d x 2_ v 2-wave gap: 



87-91) 



A, 



A^(cos k x — cos k y ). 



(11) 



We should emphasize that, even if the variational param- 
eter Ad indicates the magnitude of the d-wave singlet 
gap, a many-body state (ty d = V^d) with finite Ad does 
not necessarily mean a SC state. 92 - 93 ) To confirm the SC 
order, one needs to calculate the order parameter 13 ) or 
pair correlation functions like Pd{^) introduced later. For 
Ad = 0, $d is reduced to the Fermi sea, <&f; in this pa- 
per, we use * F (= V3>f) as the reference normal state for 
the energy gain by A^ (or Aaf)- Furthermore, we fix the 
value of t'/t in e k in the wave functions [eq. (10)] at the 
same value as that in the Hamiltonian, eq. (2), namely we 
do not explicitly introduce a band renormalization effect 
due to correlation. Because the renormalization of £k 94 ) 
brings about a serious effect in the vicinity of half fill- 
ing, 95 ) we will treat this issue in upcoming publications. 
As n goes away from half filling, the renormalization ef- 
fect rapidly becomes modest. 51 ' 88 ) 

In this paper, we consider a correlated AF-ordered 
state \I/q F (= V$af) independently of the d-wave state; 
as a one-body state, we use a Hartree-Fock solution on 
the square lattice, $af(Aaf): 96 ^ 



with 



k+Qcr' 



(12) 



(13) 



where sgn(cr) = 1 or — 1 according as o =t or I, Q is the 
AF nesting vector (tt,tt), and 



Uk(^k) = 



\ 



l-(+) 



7k 



AF 



(14) 



with 7k = — 2t(cosfca; + cos k y ). In the products of 
cq. (12), we take k-points in the Fermi sea determined 
by £k for convenience of optimization; namely, we do not 
allow for a band renormalization, similarly to 

2.3 Condition of VMC calculations 

In optimizing parameters, we employ a simple method 
that repeatedly applies one-dimensional minimization 
(Brent method) to every parameter 97 ) with the others 
fixed, because the number of parameters is small (< 6). 
In most cases, the parameters as well as energy converge 
in a few iteration loops except for special cases near Mott 
transitions. To reduce the statistical errors, we continue 
the iteration loops for more 15-20 rounds after the con- 
vergence, and average the data accumulated in this addi- 
tional process. In each loop, we use 2.5 x 10 5 (typical)-lO 6 
fixed samples; consequently, the results in optimization 
are substantially averages of several million samples, with 
typical precision of the parameters and energy being the 
order of 10~ 3 and 10~ 4 t, respectively. In the calculations 



of physical quantities using the optimized parameter sets, 
we averaged at least 2.5 x 10 5 samples. 

To check system-size dependence, we use systems of 
N s = L x L sites typically with L = 10-16, imposing 
on the periodic(x)-antiperiodic(y) boundary conditions 
to reduce level degeneracy. In principle, we aim to sat- 
isfy the closed-shell condition, because open-shell sys- 
tems bring about serious finite-size effects for small U /t's, 
especially in the SC correlation function. However, when 
we treat t'/t dependence with other parameters fixed, we 
are obliged to use systems with open shells in addition 
to those with closed shells. In this connection, as we vary 
t'/t, the manner of k-point occupation in $, especially 
in $af and $f, undergoes discontinuous changes at the 
values of t'/t specific to the system size. Consequently, 
large and irregular system-size dependence on t'/t oc- 
curs in most quantities even for large- U/t cases, as we 
will encounter in §4. At any rate, we managed not to be 
affected by special cases by checking different systems. 

3. Properties for t' — 

We discuss the results for t' = first to grasp U / 1 and 

5 dependence of various quantities. In §3.1, we consider 
the relation of the D-H binding to the Mott physics. In 
§3.2, by analyzing energies, we study the crossover of SC 
properties, competition between SC and AF, and insta- 
bility of the AF state against phase separation. In §3.3, 
we consider the U/t- dependence of correlation functions. 
In §3.4, we argue that, in the underdoped regime, 5 de- 
pendence is different between the quantities derived from 
pairing and those from charge itinerancy, and that the 
strength of SC is represented by their product. In §3.5, 
a change in the mechanism of conduction is discussed. 

3. 1 Relation of D-H binding to Mott physics 

In Fig. 2 in ref. 16, we showed that the D-H binding 
correlation is highly effective for large values of U ft and 
small doping rates. Before discussing the main results, 
here we summarize this feature for t' = 0; the feature for 
finite t'/t's is basically the same as that for t' — 0. 

Figure 1 shows the optimized D-H correlation param- 
eter fj, [eqs. (6) and (7)]. Generally, as U/t increases, \x 
increases, namely, the D-H binding becomes more tight, 
for the normal and SC states. The behavior near half 
filling is sharply distinguished according as U is smaller 
than U c (Mott critical value) or not. For U < U c , \i re- 
mains small (fi ^ 0.5), and decrease as L increases at 

6 = 0, as indicated by an arrow, whereas for U > U c , 
(i increases steadily as 5 decreases, abruptly approaches 
1 for 5 — > 0, and increases as L increases. Thus, // still 
works for the Mott physics in doped systems. For the 
AF state, the D-H factor is redundant, as discussed in 
ref. 16; the redundancy is recognized from the behavior 
almost independent of 5 in Fig. 1. 

The effect of the D-H factor is more clear in energy 
improvement. Figure 2(a) shows the difference of total 
energy per site (E) between = VqVg^cI {Eq) and 
*G = T^G^d (Eg)- The energy gain by Vq increases 
as S decreases. For U < U c , however, the gain remains 
small it half filling (see also Table I). This feature 
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Fig. 2. (Color online) Improvement of energy by D-H binding factor on simple Gutzwillcr projection for d-wave state, (a) Energy 
difference between \E"q and \Pq, as function of doping rate for five U /Vs. (b) Energy difference between and 'Pq normalized by Eq 
as function of U/t for six doping rates. The dashed lines with open symbols indicate the values obtained using >&j instead of sPq. 



Table I. The energy improvement by and over Eq of \Pq 
in percent, i.e., — Eq]/Eq x 100, is listed for three values 

of U /t and of 8. 



8 


U/t -> 


4 


8 


12 





*Q 


0.81 


48.7 


184.9 




*J 


0.08 


0.7 


5.8 


0.083 


*Q 


0.33 


17.1 


16.1 




*J 


0.06 


0.5 


0.9 


0.194 


*Q 


0.15 


2.8 


2.6 




*J 


0.08 


0.2 


0.4 



the improvement is negligible for U/t & 5. In contrast, 
the energy gain becomes sizable for U 'k, U c , where E is 
markedly improved by Vq at 5 = [Fig. 2(a)]. Although 
the energy gain by Vq decreases as 5 increases, Vq still 
brings appreciable improvement for < S ^ 0.15. This 
is roughly consistent with experiments in cuprates. 98 '") 
For large U/t's, the contribution of Vq amounts to 2/3 
of total Eq at half filling, and preserves about 7% even 
for (5-0.15. 




0.2- - 

c^o o o o — e © — s — e — g o o o~e&---o -o 

D-H binding 

0.1 ' (1.2 

S 

Fig. 1. (Color online) Optimized D-H binding parameter in the 
normal (^q), rf-wave (^q), and antiferromagnetic ( > I'q F ) states 
as function of doping rate for some U/t's. The arrows near the 
vertical axis indicate the difference of system-size dependence at 
8 = between U > U c and U < U c - The data for L = 10-16 
are plotted together. The Mott transition points arc estimated at 
U c /t ~ 7 and 9 for <]>£ and respectively. 17 ) 



The D-H binding correlation has a close relation to 
the spin exchange coupling, because a NN D-H pair 
yields a NN antiparallel spin pair with a single hopping. 
Correspondingly, the behavior of energy improvement in 
Figs. 2(a) and 2(b) becomes highly similar to that of 
the spin structure factor S'(q) at q = (n,ir) as shown 
later in Figs. 14 and 10(a), respectively. 45 ) Then, the ef- 
fective range of Vq in the U/t-S space coincides with 
the range of robust d-wave pairing [see as shown 
later in Figs. 11(a) and 5], because the d-wave pairing 
is mainly mediated by the AF spin correlation. Thus, 
the D-H binding correlation is crucial to the d-wave SC, 
and \1/q literally represents a doped Mott insulator for 
U>U C . 

Having focused on attractive correlation as intersite 
factors, here we check that intersite repulsive correlation 
factors are irrelevant near half filling. For this purpose, 
wc consider a typical NN repulsive Jastrow factor, 

^J=Pd dd Ph hh - (15) 
with pd and /?h being variational parameters, and 

Add = djd j+Tl Ahh = hjh ]+T , (16) 

where j runs over all the lattice sites, and r NN sites 
of the site j. We compare improvement in energy by 
*j = VjV G $ d (Ej) with that of ty Q over E G in Ta- 
ble I. In every case, the energy improvement by is 
much smaller than that by ^q. In Fig. 2(b), we plot 
(Eq — Ej)/Eq to compare with the case of Eq; insignif- 
icance of Vj is obvious in the present parameter regime. 
Long-range Jastrow factors are also considered irrelevant 
for S — 0, by inferring from the study at half filling. 20, 21 ) 
Henceforth, we disregard intersite repulsive factors. 

3.2 Energy gain by d-wave gap 

Now, we consider SC properties. We begin with the 
energy gain (AE) by the d-wave gap parameter A^: 

AE = E{0) - E{A° pt ), (17) 

where E(Ad) denotes the total energy for which the var- 
iational parameters other than A^ are optimized, and 
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A^ pt indicates the optimized value of A<j. Although AE 
has been often called "condensation energy", we call it 
"energy gain" or "energy difference" in this paper, be- 
cause AE does not necessarily correspond to the SC 
condensation energy measured experimentally, as will be 
discussed in §3.4. 
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Fig. 3. (Color online) (//(-dependence of energy gain AE/t for 
d-wave (solid symbols) and AF (open symbols) states for t' = 0. 
The three panels show different ranges of doping rates: (a) under- 
doped (0 < 8 < 0.1), (b) lightly underdoped to optimum-doped 
(0.1 < 6 < 0.16) and (c) over-doped (5 ~ 0.2) regimes. The AF 
state is not stabilize for S's in (c). Note that energy scales are 
different among the panels. In each panel, the data of d-wave state 
arc fitted with functions a' cxp( -3't/U) [dashed lines] in the weak- 
correlation side (U < U co ) and aexp(—U/(3t) [dash-dotted lines] 
in the strong-correlation side (U > U co ). 



Figure 3 shows the U/t dependence of AE for the d- 
wave state (AE d ) and the AF state (AE AF ) defined for 
Aaf similarly to eq. (17); the panels (a), (b) and (c) show 
the data for underdoped, optimum-doped, and over- 
doped values of 6, respectively. As discussed in ref. 16, 
AE d is negligible for U 5s 6t, but abruptly increases 
as U/t increases, and is roughly fitted as AE d /t ~ 
a'exp(-/3't/[/) 10 °) at approximately U = W{= 8t) (a', 
8': constants). AE d has a maximum at U <~ U co = 10i- 



AE d /t = aexp - — 



(18) 



12t, then decreases slowly as fitted well with a curve, 

XT 
fr, 

for 2W ^ U < AW . The estimated values of constants a 
and /3 are shown in Fig. 4 as a function of 5. From the 
form of eq. (18), we consider that the attractive pairing 
potential in this regime is proportional to t 2 /U(= <//4); 
it indicates that the low-energy physics in this regime 
is parallel to that of the corresponding i-J-type mod- 
els. 16 ' 101 ) Thus, we conclude that the properties of SC 
undergo a crossover from a BCS type to an unconven- 
tional t-J type at U ~ U co (This is a loose definition 
of U co ). As long as is SC, this behavior is qualita- 
tively independent of S, although the magnitude of AE d 
decreases with 5. 

As obvious from the form of eq. (17), AE d originates 
from a finite <i-wave gap. Actually, the behavior of AE d 
closely corresponds to that of the optimized A^ shown in 
Fig. 5; Ad abruptly starts to increase at U ~ U c , where 
the D-H binding also becomes effective, as mentioned. 

Now, we compare AE d with AE AF . As discussed in 
refs. 16 and 17, A_E AF also has a maximum at U ~ W 
(U < U co ) and has a tail for large U/t. At half filling, 
AE AF is always larger than AE d [Fig. 3 in ref. 16], al- 
though A_E AF decreases rapidly for large U/t. Almost as 
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Fig. 4. (Color online) Coefficients a (open symbols) and f3 
(closed symbols) in fitting function of AE d /t for U > U co [eq. (18)] 
as function of hole density. The circles denote the data for t' = 0. 
The data for t' 5^ are included (triangles) with the values of t'/i. 
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soon as carriers are doped, AE AF is surpassed for large 
U/t [Fig. 3(a)], and tends to vanish at a certain value 
of U/t in the optimal-doped region [Fig. 3(b)]. As 8 in- 
creases, AE AF becomes smaller than AE d (S — 0.156) 
for any value of U/t, and finally vanishes for 5 ~ 0.2 
[Fig. 3(c)]. The phase diagram constructed from these 
data [Fig. 4(a) in ref. 16] shows that the domain of d- 
wave SC tends to expand to § = and that of the AF 
state shrinks as U/W becomes large, supporting the re- 
sult of the t-J model, in which the stable state switches 
from an AF state to a d-wave state as soon as holes are 
doped. 15 ' On the other hand, the AF state seems pre- 
dominant at U ~ W. In fact, however, the AF domain 
in the phase diagram probably vanishes for 5^0 (unless 
a d- wave SC order coexists 67, 70 ' ) in accordance with the 
actual behavior of cuprates, because the AF state for a 
finite S is unstable against phase separation, as we will 
discuss shortly. 

Let us turn to the behavior of components of AE. 
In Figs. 6(a) and 6(b), the kinetic and interaction parts 
of AE d are shown as a function of U/t. For small U/t 
( 10), the energy gain is derived from the interaction 
part AEu with a loss in the kinetic part AE t . This mech- 
anism is identical with that of the BCS theory. On the 
other hand, for large U/t ( ^ 12), the SC state is stabi- 
lized by the reduction of AEt- As discussed in ref. 16, in 
this regime of U/t, the quasiparticle renormalization fac- 
tor Z in the nodal direction is enhanced by introducing 
Ad- This indicates that the SC coherence promotes the 
itinerancy of carriers, which is suppressed in the normal 
state by the strong correlation. Thus, we conclude that 
the so-called "kinetic-energy-driven SC" is realized only 
for U > U co . A FLEX calculation for the Hubbard model 
also supports this behavior, although U co /t is small and a 
sharp crossover is not obtained. 102 ' Because E t is propor- 
tional to the sum of optical conductivity <j\{uj) (exactly 
for t' — 0), 103 ' it is possible that precise optical measure- 
ments can determine whether or not the cuprates belong 
to the regime of the t-J model. 

Incidentally, similar behavior of AEu and AE t is 
found for the AF state as shown in Fig. 7, and for an 
s-wave SC state in the two-dimensional attractive Hub- 
bard model. 84 -' Accordingly, it is possibly a common ten- 
dency that nonorder-to-order transitions are induced by 
the reduction in kinetic energy in the strong-correlation 
regimes. 104 ' 

Finally, we discuss possible phase separation. This 
subject is important in relation to the density modu- 
lations or phase separation reported in the underdoped 
cuprates, 105 ' and has been theoretically pursued mainly 
for the t-J model with small J/t. 106 ' A homogeneous 
state is unstable against phase separation if the charge 
compressibility, 

1 / 2 d 2 EY 1 

K =tfX°={ n M ) ' (19) 

is negative, namely, total energy E/t is convex as a func- 
tion of electron density n. In order to check this possi- 
bility, we show, in Fig. 8, the total energy of vPq and 
4-A F as a function of S. For both U/t = 8 and 12, E/t of 



the d-wave state exhibits a concave curve, namely k > 0, 
meaning the d-wave state has intrinsic stability. In con- 
trast, the AF state exhibits a convex curve, k < 0, in- 
dicating that the AF state of partial filling gives rise to 
a phase separation of the insulating AF state with local 
filling (5 = and a SC (or normal) state with local filling 
S <: 0.2, according to the Maxwell construction. A similar 
result was reached by a recent QMC study for t' — 0. 107 ' 
This instability does not occur in electron-doped cases 
as will be discussed for t' ^ in § 4.1, and in coexisting 
states of SC and AF. 67 ' 70 ' 

3.3 U/t dependence of correlation functions 

To measure the strength of SC, a d a .2_ y 2-wave pairing 
correlation function, 

P ^)= WX £ (- 1 ) 1 ^ T ' ) x 

i r,r'— x,y 

(A£(Rj)A T / (Rj + r)) , (20) 

is appropriate for the present method. Here, x and y de- 
note the lattice vectors in the x and y directions, respec- 
tively, 5(t,t') indicates Kronecker delta, and A|(Ri) is 
the creation operator of a NN singlet pair at the site Ri , 

AUR i ) = (4 t 4 +Tl + 4 +Tt 4 l )/V2. (2i) 

If the value of -Pd(r) remains finite for |r| — > oo, & d- 
wave off-diagonal long-range order exists. As explained 
in Appendix C, we found a proper scheme to reliably 
estimate Pd(r) for |r| — > oo (P3°) from the data of finite 
L's. We confirmed that P/f° converges to zero for and 
as L increases. Henceforth, we discuss P^°. 




Fig. 9. (Color online) d-wave SC correlation function as function 
of U/t computed with d- wave states for various values of doping 
rate. Circle symbols indicate half filling (5 = 0); similarly, upward 
triangles (A) <5 ~ 0.03, squares (□) 5 ~ 0.09, downward triangles 
(V) S ~ 0.14, diamonds <5 ~ 0.20, and hexagons <5 = 0.25. The 
arrow indicates the Mott critical value at half filling. 

Figure 9 shows U / 1 dependence of P/f° calculated with 
the optimized <J/q for some 5's. Overall features for doped 
cases (5 > 0) resemble those of AE/t (Fig. 3), especially 
for the optimal and overdoped regimes (6 ^ 0.15). A large 
difference between P^° and AE/t emerges as S approach 
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Fig. 8. (Color online) Total energy as function of doping S of AF and d-wave states for (a) U/t = 8 and (b) U/t = 12. The AF state 
is reduced to a normal state for S ^ 0.18, as shown with horizontal dashed arrows. The solid lines are straight guide lines. 



rapidly as U/t increases. The extreme case is half fill- We leave an analysis of the S dependence to the next 
ing, at which vanishes for U > U c , corresponding subsection, 
to the Mott transition, as discussed in ref. 17. Thus, the 
peak value of P%°, which is a better definition of U co /t, 

rlprrpasps to trip Mnt.t. critical valnp TT~/t as A rlpnrpasps 
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For SC in the present case, the spin structure factor 



22) 



is an important quantity, because the electron pair scat- 
tering with Q = (-7T, 7r) will contribute to the d x 2_ y 2-wave 
pairing. In Fig. 10(a), U ft dependence of 5(Q) calculated 
with is plotted for six values of S. As discussed in 
ref. 17, S(C$) increases with U/t, and in particular, dis- 
continuously at U c /t at half filling in accordance with the 
first-order Mott transition. For 6 > 0, S(Q) still rapidly 
increases for U c & U & U co , although the discontinuity 
vanishes. The increase of 5(Q) for U < U co coincides 
with the behavior of P%°, but the further increase of 
S(Q) for U > U co is opposed to the decrease of P$°. This 
is because SC also depends on the itinerancy of carriers; 
for example, the quasiparticle renormalization factor Z 
in the nodal direction decreases as U/t increases [Fig. 6 
in ref. 16]. 

Similar decreasing behavior is found in the charge cor- 
relation function between NN sites, 



iVi 



TT7" ^ ( n i n i-\ 

3,T 



(23) 



as shown in Fig. 10(b). Here, r runs over the four NN 
sites of the site j. N-y represents NN charge fluctuation, 
and is related to the mobility of electrons. The decrease 
of Ni and Z, showing suppression of charge fluctuation 
by U/t, resembles the behavior of P%° for U > U c . 

The above argument implies that the strength of SC 
requires two factors, namely, the pair formation owing 
to AF spin correlation and the fluidity of pairs owing to 
charge fluctuation. We argue this topic again in the light 
of 5 dependence, which is controllable in experiments of 
the cuprates, in the next subsection. 
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Fig. 10. (Color online) (a) Spin structure factor 5(q) at q = 
(tt, it) for six electron densities. For comparison, the behavior of 
optimized d-wave gap parameter Ajj/t (exactly %A d /t + 1) is also 
shown for 8 = 0.102 (open square and dashed line), (b) Nearest- 
neighbor charge correlation function as function of U/t. The sym- 
bols for 8 and L are common between the two panels. The arrows 
indicate the Mott transition point at half filling. 
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Fig. 11. (Color online) (a) Energy gain AE/t for d-wave and 
AF states versus doping rate. The numbers near the data points 
indicate the used system sizes (L). The arrows alongside the ver- 
tical axis denote the direction of system-size dependence at <5 = 
(downward: U/t = 6; upward: U/t = 8 and 12). (b) Optimized gap 
parameters in d-wave (A^) and AF (Aaf) states as functions of 
doping rate. 



3-4 Doping-rate dependence of various quantities 

First, wc look at the energy gain of the d-wave state 
AE d [eq. (17)]. In Fig. 11(a), we show the 5 dependence 
of AE d ; AE d is the largest at half filling and monoton- 
ically decreases as S increases for any U/t. Figure 11(b) 
shows the 6 dependence of optimized d-wave gap A^, 
which is the sole parameter directly controlling the en- 
ergy scale of a singlet gap in ^q. It is natural to con- 
sider that Ad corresponds to the pseudogap. 92 ' 108 ) is 
highly similar to that of AE d as in the U/t dependence. 
Thus, we find again that is stabilized by the d-wave 
singlet formation. We will see the relation of these quan- 
tities to the spin correlation shortly. 

Next, we consider the d-wave SC correlation P/j°\ its 
S dependence is shown in Fig. 12. For U/t — 4, the mag- 
nitude of P///° is very small and has relatively large L 
and (5 dependence, indicating that firm SC is unlikely to 
appear for U/t = 4. This result is consistent with those 
of QMC calculations, 8 ' 9 ' 61 ) in which the increase in SC 
correlations as U/t increases is not found for small U /Vs. 
When U/t is just below the Mott transition point U c /t 
(~ 7), as in U/t = 6, PJ" has the maximum at half filling, 
and is basically a decreasing function of 5, except for the 
fluctuation in L (and S). This large system-size depen- 
dence suggests that the SC is still fragile. These results 
for U < U c indicates that T c does not exhibit a dome 
shape in the weakly correlated Hubbard model, although 
the approximate methods claim such behavior. 6, 109 ) 
In contrast, PJ" for U > U c displays a dome shape 
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Fig. 12. (Color online) d-wave SC correlation function versus 
doping rate for several values of U/t in ^q- The digits (10-16) de- 
note the system size L used. The stars denote A| c for U/t = 12 
calculated with eq. (24) using A<j obtained by VMC. The magni- 
tude of A| c is adjusted to be equal with the corresponding P?° . 



as a function of 8, which is caused by the vanishing of 
at half filling as a Mott insulator. The shape of P™, 
especially for U/t = 12, closely resembles that of the SC 
order parameter for the NN pairing |(cj^c^)| obtained in 
the t-J model, 13 ) and is consistent with the experimen- 
tal results of T c and condensation energy in cuprates. For 
these U /i's, SC is considered robust owing to the small 
system-size dependence. This result claims that the ef- 
fective correlation strength in cuprates be strong, i.e., 
U > U c . Thus, the high-T c cuprates are literal "doped 
Mott insulators". 1 ' 4 ' 5 ) 




Fig. 13. (Color online) (a) Quasiparticle renormalization factor 
in nodal direction versus doping rate for d-wave state, (b) Abso- 
lute values of and E^ are shown as functions of 8 by solid 
symbols in a strongly correlated case (U/t = 32). For comparison, 
we add the doublon density (open triangle, right axis), a guide line 
oc C(8) = 25/(1 + 8) (dashed line), and the energy of free electrons 
as a function of electron density n (open circle). 

For U > U c , increases almost linearly to 8 in the 
underdoped regime, in contrast with the quantities like 
Ad/t and AE d /t, which are monotonically decreasing 
functions of 8. As mentioned, this is because the strength 
of SC depends on two factors, singlet-pair formation ( Ad) 
and quantities related to charge transport like carrier 
density and mobility of carriers. The latter quantities are 
bound to increase as 8 increases for U > U c , released from 



As an example, in Fig. 13(a), we show the quasiparti- 
cle renormalization factor Z estimated from the jumps 
in rt(k) in the nodal direction. Since Z roughly repre- 
sents the inverse effective mass, Z is zero at half filling 
and monotonically increases as 8 increases for U > U c , 
in accordance with the result for the t-J model. 108 ) N\ 
[eq. (23)] and the conductive part of E t discussed in §3.5 
[\Eh\ in Fig. 13(b)] are also increasing functions of 8. 

To consider 8 dependence of SC strength, it is useful 
to mention a simple analytic calculation. In a pioneering 
study using a Gutzwiller-type approximation for the t-J 
model, 92 ) a relation 

A sc = ^A,. (24) 

was derived; here the energy scale of SC (Asc) originates 
solely in A^ (d-wave singlet gap) , but modified by the fac- 
tor of carrier density. As an example, we estimate Ag C 
for U/t = 12 using A d calculated by VMC, and the re- 
sult is plotted in Fig. 12 by open stars. Ag C qualitatively 
agrees with P^° for U/t = 12. A similar dome shape is 
also obtained from ZAd, instead of eq. (24). Note that 
in a slave-boson- mean- field theory, 4 ' no ) T c in the under- 
doped regime is determined by the Bose-condensation 
temperature Tb (oc 8) of holons, which represent the 
charge part of carriers. 




S 



Fig. 14. (Color online) Spin structure factor at q = (w, vr) in 
d-wave state as function of doping rate for five U/t's. The inset 
represents S(q) at U/t = 12 for various doping rates along the 
path (0,0) -> (it, 0) -> (7T, 7r) -> (0,0). The dashed line shows the 
value for U = and 5 = 0. 

We turn to the spin structure factor 5(q). In the in- 
set of Fig. 14, S(q) of ^ d Q at U/t = 12 is depicted for 
various 8's. A sharp peak at q = Q near half filling con- 
firms a predominant AF spin correlation, and S(q) pre- 
serves the maximum at Q for 8 0.2. The main panel 
of Fig. 14 shows the 8 dependence of S(Q) for five val- 
ues of U/t; S(Q) is a decreasing function of 8 for any 
correlation strength. The small system size dependence 
indicates a short-ranged nature of 5(Q). We here repeat 
the close connection between the spin correlation and the 
D-H binding correlation (Fig. 2); 8 dependences of AE d , 
Ad and 5(Q) closely resemble to one another. This sug- 
gests a strong correlation between A^ and AF spin cor- 
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by the formation of d-wave singlet pairs through the AF 
spin correlation. 
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Fig. 15. (Color online) Small-|q| behavior of (a) spin and (b) 
charge density structure factors in (0,0)-(7r,0) direction for U/t = 
12 and some doping rates. In (b), JV(|q|) oc |qp with 7 < 1 for 
q — > seems to hold for 8 > 0. The system used and symbols of S 
are common to the two panels. 



Now, we touch on how spin and charge density gaps 
evolve when 5 is introduced. It was revealed in ref. 17 
for half filling that "Jq has a finite spin (SC) gap for 
any positive U /t, but is gapless for U < U c and gapped 
for U > U c in the charge sector. It is known within the 
single-mode approximation 111 ' that a charge density gap 
opens [closes] if the charge density structure factor, 

N W = W E e iq ' (R '- R3) M - n\ (25) 

behaves as iV(q) oc |q| 2 [oc |q|] for |q| — > 0. For the 
spin sector, a similar treatment with S*(q) is available. 
In Fig. 15, we show small- |q| behaviors of S(q) and 
N(q) at U = 12t (> U c ) for various doping rates. A 
SC (singlet) gap survives up to the overdoped regime 
{5 > 0.15), although the quadratic feature of <S(q) be- 
comes less distinctive as S increases. On the other hand, 
in -/V(q), the quadratic behavior immediately disappears 
on being doped, namely, the charge density gap vanishes 
for S > 0, because |q| 2 /A(q) seems to vanish for |q| — > 0. 
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Fig. 16. (Color online) Schematic <5 dependence of two kinds of 
gap scales deduced from present calculations for U > U c . 



present calculations (Fig. 16). The gap scales represent- 
ing the features of SC for U > U c are classified into two 
kinds: (i) Quantities exclusively related to the singlet for- 
mation, which is derived from the superexchange interac- 
tion, monotonically decrease by doping, which weakens 
the AF correlation. They are symbolically indicated by 
Ad in Fig. 16. (ii) Quantities directly related to SC, which 
is affected by both the singlet gap and charge trans- 
portability, increase as 6 increases in the underdoped 
regime, as shown by Asc- It is natural to substitute T c 
(T pair ) for A sc (A d ) in Fig. 16. 

Finally, we touch on the relation between AE d in 
eq. (17) and the SC condensation energy experimentally 
observed (Ai? con( j). In the BCS theory, AE d coincides 
with AE con d at T = 0. However, this feature is not 
applicable to the present case, because only a part of 
Ad contributes to Asc and the residual part remains 
as an incoherent singlet gap. Actually, A_B con d's esti- 
mated from the specific heat measurement 112 ' 113 ' exhibit 
a dome shape as a function of d as in T c , and is utterly 
different from the monotonic behavior of AE d shown in 
Fig. 11(b). Correctly, the above incoherent part should 
be cancelled with the corresponding part of a proper nor- 
mal state; improvement of the normal state is necessary. 
Anyway, we should be prudent in comparing AE d with 

A-Econd- 

3.5 Mechanism of conduction 

Here, we show that the mechanism of conduction un- 
dergoes a marked change through the crossover. To this 
end, it is useful to analyze the kinetic energy Et into two 
components, Et = Ed + Eh, according as the hopping 
process varies (Ed) or does not vary (£"h) the number of 
doublons, 114 ' as shown in Fig. 17(a). Ed is derived from 
the hopping processes of creating and destroying D-H 
pairs, and corresponds to the J term in the t-J model. E^ 
comes from the direct hopping of holons and doublons, 
which corresponds to the t term in the t-J model. 115) In 
Fig. 17(b), we plot E^ and E^ as a function of U/t for 
several 5's. The behavior of the two components changes 
definitely at approximately U co /t. 

We begin with the strongly correlated regime (U > 
U co ). At half filling, E^ substantially vanishes in the Mott 
insulating regime (U > U c ); Ed remains finite and be- 
haves proportionally to — At 2 /U (= — J). Since the state 
is insulating, Ed here does not contribute to current, 
namely, local processes of creating and annihilating D-H 
pairs are only repeated. When carriers are doped, this 
behavior of Ed (oc — J) is basically unchanged; its mag- 
nitude decreases slowly as S increases, accurately corre- 
sponding to the decrease of doublon number, as shown 
in Fig. 13(b). Thus, the local D-H processes at 6 = re- 
main intact for S > 0, meaning that Ed is not involved in 
conduction. On the other hand, E^ becomes finite and 
is still almost constant for U > U co with the magni- 
tude linear in S or oc 25/(1 + 5) in eq. (24), as depicted 
in Fig. 13(b). This indicates that independent motion 
of doped holons is realized for U > U co , although their 
mass is somewhat heavier than that of the free electrons. 
It follows that two kinds of holons play entirely differ- 
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Fig. 17. (Color online) (a) Hopping processes which contribute 
to the two components of kinetic energy and E^ are schemat- 
ically shown, (b) U/t dependences of E d and E^ are shown for 
six doping rates. The two dash-dotted lines in a small-f//t [large- 



. . o»ooo 
(a) ooooo 
tt tt oo»oo 
f < f n ,ooooo 
• oooo 
ootoo 
ooooo 
oooot 
otooo 
ooooo 
oooo* 
o»ooo 
ooo»o 



o»oto«ooo«oo 

oooooooooooo 
•ooo*c*«oooo 

0000*0000*00 
0000««0»0000 

o»oo#oootooo 
ooo»o«ooooo» 

000000000*00 

ooo*oooo»oo« 
oootoiooooo* 

oooooooootoo 
ooo»oooo»ooi 

ooooooiooooo 



... ooooooooooooooooo 

(b) 0000*0000000000*0 

ooooooooooooooooo 
t /> t / coOOOOOOO o'«vo o o o o o o o 
o o o o o o o Oo,o oooo o*o 
00*00000800000000 
ooooooooooooooooo 
ooooooooooooooooo 
o o o o b • » o o o o o 0*0 o o 

OOOOOTJOOQOOOOOOOO 
OOOOOOO Q'OVD 0000000 
0000000 0»,0 OOOOOOO 
000*000000000*000 



oooo»oo 

OOOOOOO 

oo»o§*o 

• 90®«00 

009««00 

o»oooo» 

OOOOOOO 

ooo«ooo 

OOOOOOO 
OOOOOOO 

ooo»ooo 

OOOOOOO 

oo»oo»o 

OOOOOOO 
OO »DOOO 

OOfiOOOO 
OOOOOOO 

000000* 

OOOOOOO 
OOOOOOO 
OOOOOOO 
OOOOOOO 

000*000 

OOOOOOO 

o o o'« d» o 

0000C00 



• doublon 

O holon 
O spinon 




Fig. 18. (Color online) Mechanism of conduction is schemati- 
cally compared between (a) the weakly correlated [Fermi liquid or 
BCS-type SC] regime and (b) the strongly correlated [doped Mott 
insulator] regime. Small arrows represent the motion of particles if 
a field (E) is applied; in (a) arrows are drawn only for a small part 
of particles for clarity. The dashed ellipses in (b) indicate the local 
creation and annihilation processes of bound D-H pairs, which do 
not contribute to current. 



U/t] regime are guides as ~ (U/t) 
[~ -t/U = - J '/(At) in E d for 8 = and 0.0833]. 



const, for 8 



0.1389 



devote themselves to forming local singlet pairs, whereas 
the holons introduced by doping work as current carri- 
ers, as schematically sketched in Fig. 18(b). Thus, the 
number of doped holes is equal to the number of carri- 
ers, which will make a small quasi Fermi surface like the 
Fermi arc or pocket. This feature is nothing but that of 
the t-J model, 5 ) and consistent with various experiments 
for the cuprates, 116 ) and supports the very low superfluid 
densities. 117 * ) The residual (background) spinons (singly 
occupied sites) remain localized unless the doped holons 
collide with them; this fact is possibly related to a recent 
neutron scattering for Bi2212, 118 ) which argued that the 
source of magnetic response in doped cuprates is local- 
ized spins. 

In contrast, for U < U co , both E d and Eh behave as 

E d ~ c d U 2 /t + E d (0), 

E h ~ c h U 2 /t + E h (0), 

with Cd and c n being constants, as shown in Fig. 17(b). 
This result indicates that every hopping process gives a 
contribution to E t in the same way. Since the D-H bind- 
ing is ineffective in this regime, holons cannot be clas- 
sified into the two kinds, and all electrons contribute to 
conduction, as in Fig. 18(a), resulting in a large (ordi- 
nary) Fermi surface with the carrier number of N . Cor- 
respondingly, the number of doped holes discords with 
the carrier number, contrary to the experimental results 
in cuprates. 

We have also obtained results similar to Fig. 17(b) for 
the normal and AF states, although we do not show here. 
The mechanism of conduction in doped Mott insulators 
entirely differs from that in weakly correlated systems. 



4. Effect of diagonal transfer t' 

In this section, we study the effect of next-nearest- 
neighbor hopping {t'), which is the principal term to 
characterize the individuality of each high-T c cupratc. 
In §4.1, we consider the stability of against ^q F in 
introducing t'and possible phase separation near half fill- 
ing. In §4.2, we discuss the origin of stability of d-wave 
and AF states on the basis of the effect of t'. In §4.3, we 
study the effect of t'/t on P$°, whereby we discuss the 
mechanism to enhance P5°. 



4-1 Energy gain by d-wave gap 
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Fig. 19. (Color online) Energy gain for d-wave (solid and half- 
solid symbols) and AF (open and double-line symbols) states as 
function of t'/t for optimally-doped densities. For comparison, data 
for some values of U/t and two values of L are plotted. 
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because there is a symmetry between t' and —t', as dis- 
cussed in ref. 17 and Appendix A. On doping carriers, 
however, this symmetry is broken. Figure 19 shows AE d 
and A_E AF for some lift's as a function of t'/t at roughly 
optimally doped densities S = 0.12 and 0.139. For the 
hole-doped cases (t'/t < 0), AE d tends to be enhanced, 
and AE AF is suppressed. For U/t > 8, AE d has a broad 
maximum in the range —0.3 t'/t — 0.1, and tends 
to vanish relatively rapidly for t' ft <> — 0.45. In contrast, 
AE AF vanishes except for a small peak at t'/t ~ —0.15. 
This peak is caused by the overlap of the Fermi sur- 
face with flat-band parts near (%, 0), as we will discuss 
later for \l/g. For the electron-doped cases (t'/t > 0), the 
situation is opposite. AE d gradually and monotonically 
decreases as t'/t increases, whereas AE AF monotonically 
increases until sudden breakdown at a large t'/t. 

The U/t dependences (not shown) of AE d and AE AF 
are qualitatively similar to those for t' = in Fig. 3, 
although the overall amplitude depends on t'/t as in 
Fig. 19. 
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Fig. 20. (Color online) (a) Kinetic and (b) interaction parts of 
energy gain AE d for various values of t'/t as function of U/t. 

We touch on the components of the energy gain AE 
for finite t'/t. For small \t'/t\, the sum of cr 1 (w) becomes 
approximately connected to the kinetic energy -Ekin- In 
Figs. 20(a) and 20(b), we plot the kinetic (t and t' terms) 
and interaction parts of AE, respectively, as functions of 
U /t. The behavior of the two components are qualita- 
tively identical to that for t' = (Fig. 6): SC is induced 
by the reduction in interaction (kinetic) energy for small 
(large) U/t's. Here, the doping rate is fixed at 5 — 0.12, 
but the tendency is unchanged for other 5's. 

In §3.2, we showed that ^q F is unstable against phase 
separation for t' = (Fig. 8). Here, we discuss the ten- 
dency toward the phase separation for t' ^ 0. In Fig. 21, 
we plot the total energy of the AF and d-wave states as 
functions of 8 for hole- and electron-doped cases. For ^q, 
E d /t is concave (k > 0) for both values of t'/t, indicating 
that the d-wave state is still stable if t' is added, regard- 
less of the kind of doped carriers. For the AF state, it 
becomes slightly concave for t'/t — 0.2 (electron-doped 
case) [Figs. 21(b)], indicating that the commensurate AF 
state is intrinsically stable for the electron doping. On the 
other hand, E AF /t is still convex (k < 0) for t'/t = -0.10 
(hole-doped cases) in the whole range of finite optimized 



a phase separation occurs in the underdoped regime as 
in t' = 0. Such behavior is almost independent of the 
value of U/t as far as U/t is sufficiently large; the same 
behavior has been found for the t-J model. 119 ) 

The reason of this instability of the AF state seems 
trivial; doping rapidly deteriorates the nesting of the 
Fermi surface for t'/t < 0, in contrast with the case 
of t'/t > 0. Although we have to be careful to com- 
pare directly the present result with experiments, be- 
cause the band renormalization effect is not explicitly 
introduced in *q and \&q F , there are some sugges- 
tive experiments: A neutron scattering experiment in 
very lightly doped La 2 - a; Sr :r Cu04 with x < 0.02 120 ) 
found a phase separation to a commensurate AF ordered 
phase and a spin-glass phase. On the other hand, such a 
phase separation does not appear in the electron doped 
Pri_ a; LaCe x CuO4. 120 ' 121) 

4-2 Origin of t' -dependence regarding stability of d- 
wave and AF states 

The t'/t dependence of AE d and AE AF in the opti- 
mally doped regime (Fig. 19) is basically unchanging if 
S varies, as shown in Fig. 22 for U/t = 12. For t'/t = 0.2 
(solid symbols), we have AE AF > AE d for any 5; for 
t'/t = (half solid symbols) and —0.1 (double-line sym- 
bols), AE AF and AE d interchange at 5 ~ 0.89 and 
0.91, respectively, whereas for t'/t = —0.25 (open sym- 
bols), we always have AE AF < AE d . Thus, the range 
of 6 where vPaf (^d) is predominant rapidly shrinks as 
\t' /t\ increases on the hole- (electron-)doped side. In- 
cidentally, we find AE AF remains small near half fill- 
ing for t'/t = 0.2 and —0.25. This is probably because 
the renormalization of the Fermi surface is not explic- 
itly introduced in ^q F ; this effect becomes serious for 
5 — » to retrieve the nesting condition. 122,123 ) We ex- 
pect AE AF > AE d for S ~ even for large \t'/t\ in 
improved calculations. 

On the basis of Figs. 19 and 22, let us think about the 
stability of d-wave and AF states with respect to t'/t. 
As \t' /t\ increases, the curvature of the Fermi surface in 
$f becomes more concave in the nodal direction (0,0)- 
(tt, 7r) as shown in the inset of Fig. 22. Furthermore, it is 
expected that the scattering with q = Q is activated in 
the area of doped Mott insulators, as discussed in §3.3 
and §3.4. 

First, we discuss the origin of stability of x J'q F . In the 
electron-doped cases, the nesting does not deteriorate 
rapidly when S increases away from half filling, namely, 
the Fermi surface continues to largely overlap with the 
magnetic Brillouin zone boundary as a whole, as shown 
in the inset of Fig. 22. The effect of large density of state 
near (it, 0) [see Fig. 23] is subsidiary for *q F - On the 
other hand, the Neel order disappears rather suddenly, 
when t'/t becomes excessively large or holes are doped, 
because the nesting condition becomes insufficiently sat- 
isfied. Since the band renormalization effects becomes 
important for the AF order for U > VF, 88,123 ) we will 
reconsider it in future publications. 

In contrast, conceivable reasons why the d-wave favors 
finite negative values of t'/t are as follows: (i) Because the 
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Fig. 21. (Color online) Total energy of d-wave (^q) and AF (^E'q F ) states as function of doping rate. In (a) and (b), a hole-doped 
{t'/t < 0) and an electron-doped (t'/t > 0) cases are plotted, respectively. The data for several system sizes are fitted by the method of 
least squares, and shown with dashed (AF) and dash-dotted (d-wave) lines. The gray solid lines are straight guide lines for the AF case. 
In the region indicated by 'normal', the optimized AF state is reduced to the normal state, namely, Aap = 0. In (b), we omit the data 
of ^q F in the very vicinity of half filling owing to an abrupt phase change in ■ Similar data for t' = are given in Fig. 8. 



(i-wave gap is minimum in the nodal directions and max- 
imum in the antinodal (k x and k y ) directions, the energy 
gain is large when the electrons near (0, tt) and (tt, 0) are 
occupied, (ii) k-dependent density of state, l/|Ve(k)|, is 
minimum in the nodal direction including (tt/2, tt/2), but 
becomes large near the antinodal points (tt, 0) and (0, tt) 
owing to the band flatness. Actually, l/|Ve(k)| becomes 
large for t'/t = —0.3 as q approaches (tt, 0) from (0,0), 
as in the dashed lines in Fig. 23, which shows l/|Ve(k)| 
in the noninteracting case, (iii) The vector Q connects 
two k points near two nearest antinodal points with mu- 
tually inverse signs of Ak- Thus, can take advantage 
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Fig. 22. (Color online) Comparison of energy gain between d- 
wave (triangles) and AF (circles) states as function of doping rate 
for four t'/t's and U/t = 12. The arrows indicate the crossing points 
of AB AF and AE d for t'/t = and -0.1. Systems of L = 10-16 
are used. The inset shows the Fermi surface in the electron picture 
for n = 0.85 (favored by d-wave) and 1.10 (favored by AF) both 
for t'/t = —0.3. The arrows are the Q vectors connecting the hot 
spots [intersections of the Fermi surface and magnetic Brillouin 
zone boundary indicated by the gray bold line]. The gray dotted 
lines indicate the nodes of the d-wave state, and the gray circles 
the van Hove singularity points near the Fermi surface. 



of d-wave pair scattering of Q near the antinodal points, 
whereas the scattering around (tt/2, tt/2) seems of little 
importance. We will pursue this topic in §4.3 and §5. 

This feature is corroborated by the behavior of the 
momentum distribution function, 

n(k) = i^(cLc kff >. (26) 

Figure 23 compares n(k) among the optimized (i-wave, 
AF and normal states for an underdoped density. n(k) 
of the normal state exhibits clear discontinuities (Fermi 
surface) near (tt,0) and (tt/2, tt/2). Near (tt,0), n(k)'s 
of both (i-wave and AF states are considerably modi- 
fied, meaning scattering actively takes place there. Near 
(tt/2, tt/2), n(k) of the AF state is again noticeably 
changed from that of the normal state depending on t'/t, 
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Fig. 23. (Color online) The momentum distribution functions 
are compared among the d-wave (open symbols) , AF (half closed) 
and normal (closed) states for a couple of t'/i's, S = 0.094 and 
L = 16. Firm orders exist for the d-wave and AF states. The gray 
horizontal line shows n/2. The gray dashed curves represent the 
relative magnitude of noninteracting k-dependent density of states 
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but n(k) of the c?-wave state changes only slightly, and is 
almost independent of t' jt. Thus, for ^q, the scattering 
in the nodal direction is ineffective, but the electronic 
state near the antinodal points should be noticed. 

4-3 d-wave pairing correlation 

Now, we turn to the d-wave SC correlation function 
in ^q. In Fig. 24, U/t dependence of Pf , defined by 
eq. (20) and in Appendix C, is shown for four finite t' /Vs. 
In Fig. 25, S dependence of P£° is shown for four U/Vs. 
Although the overall behavior is similar to the case of 
t' = (Figs. 9 and 12), there are noteworthy differences 
owing to the effect of V . (i) The magnitude of P£° is 
somewhat enhanced for t'/t = —0.1 and —0.25, but is 
suppressed for t'/t = 0.2, in accordance with the behav- 
ior of AE d discussed in §4.1. (ii) For U > U c , the doping 
rate giving the maximum P^° for a fixed value of U/t 
shifts to the overdoped side as t'/t decreases, (iii) For 
t'/t = —0.1, P//° starts to increase at a smaller value 
of U/t. (iv) Although we do not show detailed data in 
Fig. 24, for t'/t <, - 0.4 and 5 ^ 0.2, a state with large 
PJ° is competitive with a state with tiny P£° (A^/t ~ 0), 
suggesting that the SC state collapses at approximately 
these parameters. 

First, let us consider the point (i) and (ii) in the 
strongly correlated regime. As in Fig. 25(d), for U (= 12t) 
sufficiently larger than U c , P//° forms a well-proportioned 
dome shape as a function of S. This S dependence of P£° 
is consistent with that of T c and condensation energy 
observed in cuprates. We define <5 max as the value of 6 
which marks the largest P^° for given U/t and t'/t. For 
S < <5max, P//° is proportional to 6, and is almost inde- 
pendent of the values of t'/t and L. This means that SC 
in the undcrdoped regime is steady, and the magnitude of 
P///° depends only on U/t, and not on minute band struc- 
ture. On the other hand, for 5 > <5 max for the hole-doped 
cases (t'/t < 0), P//° irregularly depends on the values 
of t'/t and L, and is scattered to some extent. 124 ) Thus, 
the strength of SC is sensitive to the shape of Fermi 
surface and fragile in the overdoped regime. However, 
P^° becomes larger than that for t' = 0, probably be- 
cause the t' term bends the Fermi surface so that it may 
pass by the antinodal points. Similar behavior has been 
known for the t-J model. 51,88 ' For the electron-doped 
case, <5 max is smaller than those for the hole-doped cases, 
but the fluctuation with respect to S and L is small for 
S > Smax, so that SC is steady. In the transitional regime 
[Fig. 25(c) for U c < U & U co ], although P//j° manages 
to preserve the dome shape, S max is lowered to have a 
wide irregular range for t'/t < 0. For the electron-doped 
case (t'/t = 0.2), P^° is suppressed as compared with 
t'/t = 0, and almost vanishes in the overdoped regime. 
Some similar features were pointed out in a recent study 
by Tocchio et a/. 125 ) 

Next, we look at the case of intermediate correlation 
strength [U = 6t, Figs. 25(b)], which is slightly smaller 
than the Mott critical value U c . As discussed in ref. 17, 
Pj° for t' = exhibits a sharp peak immediately below 
U c (~ 7t) at half filling; such behavior is unchanged for 
t' 7^ (black circle symbols in Fig. 24). Because the value 



largest at half filling, and tends to decrease as S increase. 
Especially for t'/t = —0.4 and 0.2, P/j° almost vanishes 
for 5 > 0.05. Similarly to the case of t' = 0, this behavior 
of P£° is contradictory to the dome shape of T c and con- 
densation energy observed in the cuprates, in contrast 
with the cases of U > U c . 

Now we move to the point (iii) regarding small values 
of U/t. As discussed in §3.3 for t' = 0, PJ° is negligibly 
small for U/t as small as 4. As found in Fig. 24, this fea- 
ture basically does not alter even if t' is added, except 
for t'/t = —0.1 [Fig. 24(a)]. This exceptional enhance- 
ment of P/j° is evident in Fig. 25(a), where appreciable 
magnitude appears only for t'/t = —0.1 in the six values 
of t'/t. We argue in the following that this exceptional 
enhancement of P/j° at t'/t ~ —0.1 for a small U/t is 
useful to infer the origin of SC. 

In Fig. 26, we show P^° for U/t = 4 with star sym- 
bols as a function of t'/t for four <5's. For an under- 
doped density, S = 0.0833 [Fig. 26(b)], the meaning- 
ful magnitude of PJ° is limited to a narrow range of 
—0.16 t'/t ^ —0.075. Similarly it is also limited to the 
range of —0.275 & t' /t — 0.075 for an overdoped den- 
sity, S — 0.1944 [Fig. 26(d)]. These results are consistent 
with that of a recent QMC study, 61 ) which concluded 
that the SC susceptibility vanishes for t' = but remains 
finite for t'/t = —0.2 for U/t < 5 in the optimally doped 
regime. In Fig. 26, we simultaneously plot the momentum 
distribution function n(k) calculated with \&q and $>q at 
a couple of available k points near the antinodal point 
tt(1, 0), i.e., X = tt(1, i) and X' = tt(1 - f , ±). As t'/t 
is increased, n(X) or n(X') of iPq discontinuously drops 
from near unity to near zero, when the Fermi surface 
passes through the X or X' point, as indicated by arrows 
on the upper axes. Note that in each panel, the positions 
of the peak and shoulders of P^° precisely coincide with 
those at which the Fermi surface overlaps with antinodal 
k points. As in Fig. 23, (±7r, 0) and (0, ±7r) are van Hove 
singularity points for \t' /t\ < 0.5, and connected one an- 
other by the AF vectors Q = (ir,±ir), (±tt,tt). It fol- 
lows that the pair scattering with Q and SC thereby are 
enhanced at the values of t'/t indicated by the arrows. 
This SC mechanism is consistent with the c?-wave BCS 
theory. 126 ) 

However, this mechanism for small U/t's has difficul- 
ties in explaining features of the high-T c cuprates. First, 
the effective values of t'/t of the cuprates show a rela- 
tively wide range; single-layer La systems and double- 
layer Y and Bi systems have t'/t ~ —0.1 and —0.3, re- 
spectively. 46 ~ 50 ) However, the above SC mechanism sen- 
sitively depends on t'/t or the band structure, and SC 
correlation is enhanced only in the limited range. It is 
unlikely that robust SC occurs at t'/t = —0.3 in the 
whole relevant range of 5. Second, it is probable that 
the magnitude of P%° discussed above is, as a matter of 
fact, much weaker. As shown in Fig. 25(a), P£° has a 
large and relatively regular system-size dependence, and 
consequently seems to become negligible for L — > oo. 

For comparison, we make the same analysis for a 
strongly correlated case, U/t — 12 (Fig. 27). In the un- 
dcrdoped regime [panels (a) and (b)], P%° has apprecia- 
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Fig. 24. (Color online) d-wave pairing correlation function as function of U/t for various values of doping rate 8 (rounded off to 
two decimal places in legends). The value of t' /t is different among the four panels: (a)-(c) hole-doped and (d) electron-doped cases. 
Corresponding results for t' /t = are given in Fig. 9. The arrow in each panel indicates the Mott transition point at half filling (circles). 




Fig. 25. (Color online) d-wave pairing correlation function as function of doping rate 8 for various values of t' /t. The value of U/t is 
different among the four panels (a)-(d). The symbols indicating the values of t'/t are common to the four panels. Small digits beside 
some data points indicate the used lattice size L (= 10-16). 



pcndence is very weak, even between electron-doped and doping rate [Fig. 27(c)], the electron- hole asymmetry 
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Fig. 26. (Color online) Correspondence between d-wave SC correlation function and jump in momentum distribution functions of vT/q 
and at two k points (k = X and X') near (tt, 0), as t'/t varies, for a weak correlation (U/t = 4). The cases of four doping rates are 
shown for L = 12. The arrows on upper axes indicate the values of t'/t at which the Fermi surfaces of pass X and X'. 
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Fig. 27. (Color online) The same quantities as those in Fig. 26 are depicted in the same scales for a strongly correlation (U/t = 12). 
The band renormalizaion by correlations is not considered here, but its effect seems irrelevant except near half filling. 51, 88, 94 ' 123 ) 



in the electron-doped regime t'/t > 0, but P^° is steady 
and large throughout the range of t'/t < 0. This tendency 
is preserved for an overdoped regime [Fig. 27(d)], where 
P%° is still large in the range of the hole-doped cuprates 
(-0.4 t'/t ^ - 0.05). On the other hand, the Fermi 
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the weakly correlated cases as indicated by arrows. For 
the underdoped and optimally-doped cases, P%° does not 
exhibit any special behavior at the positions of the ar- 
rows. Only in the overdoped case, P%° displays a slight 
tendency of forming a shoulder near an arrow. Anyway, 
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most independent of the locus of the Fermi surface of the 
underlying normal state, namely, the bare band struc- 
ture, suggesting that P//° is determined by a mechanism 
qualitatively different from conventional BCS-type theo- 
ries, which start from the instability of the Fermi surface 
of the normal state. We return to this subject in §5. 

To summarize the t'/t dependence of PJ* 5 , we show 
P//° for three different regimes of 5 in Fig. 28. For a 
weak correlation (U/t — 4), the area of enhanced P^° is 
limited (-0.25 <> t'/t <, - 0.1) regardless of 6, and the 
magnitude tends to vanish as L increases. For strongly 
correlated cases (U/t > 8), P£° is enhanced by moderate 
negative values of t'/t, when <5 is in the overdoped and 
optimally doped regime. On the other hand, near half 
filling, the magnitude of P^° is steady, and tends to be 
independent of t'/t, as U /t increases. 

Finally, we point out that P%° becomes a weakly in- 
creasing function of t'/t near half filling and for very large 
U/t, as shown in Fig. 28(d). Such unexpected behavior 
was discovered for very slightly doped t-J models, using 
density matrix renormalization group 127 ' and exact diag- 
onalization, 128 ) and seemed incompatible with the prop- 
erty of cuprates. Afterward, a VMC study 51 - 1 for the t-J 
model revealed that this feature is restricted to the very 
vicinity of half filling. The present result supports this 
finding, and add a requirement that the interaction be 
considerably strong (e.g. U/t = 16). 

5. Antinodal Electrons and Superconductivity 

In §4.2, we argued that the electrons near the antinodal 
points are crucial for the d-wave SC. In this section, we 
actually reveal a close relationship between the behavior 
of momentum distribution function n(k) near k = (-7T, 0) 
and the SC correlation function P^° in the strongly cor- 
related regime. 

First, let us recall again the t'/t dependence of P^° for 
sufficiently large U /t's to compare it with that of n(k). As 
discussed in Ref. 17, properties of in the insulating 
phase (U > U c at half filling) are almost independent 
of t'/t. For small S, the t'/t dependence of P|° is still 
weak as seen for U/t — 12 and 16 in Fig. 28(c), but for 
S <i 0.1, P/j° comes to largely depend on t'/t, i.e., P/j° 
tends to increase for t'/t < and decrease for t'/t > 
as in Figs. 28(b) and 28(a). 

With this behavior of P^° in mind, we look at the 
t'/t dependence of n(k) [eq. (26)]. In Fig. 29, n(k) for 
six S's is plotted for U/t = 12. In each panel, data for 
various values of t'/t are displayed together. At half fill- 
ing [Fig. 29(a)], n(k) is almost independent of t'/t as 
mentioned, and does not have a discontinuity at any 
k, because is Mott insulating (U c /t ~ 6.5-7.2) with 
a charge density gap. On doping carriers, a discontinu- 
ity (Fermi surface) appears in the lattice-diagonal direc- 
tion near (ir/2,%/2) [Figs. 29(b)- 29(f)], 108 > because ^ d Q 
turns to SC with nodes of Ak for k x = ±k y (and fi- 
nally metallic). Meanwhile, it is in the antinodal area 
near (tt, 0) where the behavior of n(k) markedly changes 
as t'/t varies; the degree of change culminates around 
the optimum doping rate (5 ~ 0.15). Compared with the 
antinodal area, the change in n(k) by t'/t is insignifi- 



close correspondence between SC correlation (Pj°) and 
the electronic structure in the antinodal area. 

In the overdoped regime (S ^ 0.2), n(k) exhibits 
a Fermi-liquid-like discontinuity (Fermi surface) in the 
antinodal area for some t' /t's, i.e., on the segment (0, 0)- 
(ir, 0) in the electron-doped cases (t'/t > 0), and on 
(7r,0)-(7r,7r) for t'/t £ - 0.4 [Figs. 29(e) and 29(f)]. We 
find in Fig. 28(a) that when such a discontinuity appears, 
robust SC does not occur. Inversely, when n(k) changes 
continuously and relatively slowly in the antinodal area 
[Figs. 29(b)-(d)], P%° has a large and steady magnitude 
[Figs. 28(b) and (c)]. Such behavior appears in a broad 
range of t'/t particularly in the underdoped regime. We 
will discuss this problem shortly. 



Antinodal Nodal Antinodal Nodal 




(0,0) (71,0) (TC,7C) (jt/2,jt/2) (0,0) (71,0) (7i/2,7i/2) 

k k 

Fig. 30. (Color online) Elements of bare band dispersion rela- 
tions along paths (0, 0) — > (tt, 0) — > (tt, tt) — > (0, 0) (left panel) 
and (7r,0) — > (7r/2, 7r/2) (right panel): e/t = — 2(cosfe a; +cosfe s ), 
e'/t' = -Acosk x cosk v , e"/t" = — 2(cos2fc;c + cos 2k y ). t" indi- 
cates the hopping integral to the third-neighbor sites (±2, 0) and 
(0,±2), which is disregarded in this paper. The nodal (7r/2,7r/2) 
and antinodal (71-, 0) points are marked by vertical dashed lines. 



The above feature of rc(k) is more or less related to the 
band structure. In Fig. 30, the elements of the bare band 
dispersion relation £k [= e + e'(+e")) is depicted. The 
dispersion of the pure square lattice e (f = t" = 0) has 
a well-known flat part at e = around (tt, 0) (stationary 
in both k x and k y directions), where the Fermi surface 
passes for <5 ~ 0. The dispersion due to t' (e') is also sta- 
tionary but maximum at (71", 0). Consequently, and the 
locus of Fermi surface around (tt, 0) sensitively depend on 
t'/t. In contrast, at the e = point in the nodal direction 
(7r/2,7r/2), e' is zero and stationary, so that the change 
of £k by t'/t is slight near (tt/2,tt/2). Incidentally, the 
third-neighbor hopping t" , which is omitted here but is 
often used as t" = —t'/2 in the literature, works similarly 
to that oit' near (tt,0) because e"/t" e'/t'. 129 *) 

To relate the above feature of n(k) to P/j° quantita- 
tively, we need a quantity that appropriately expresses 
this feature. The magnitude of n(k) near (tt, 0) itself ap- 
parently does not represent the strength of SC. In fact, 
we found that the slope of n(k), 
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Fig. 28. (Color online) d-wave pairing correlation function as function oft' /t for several values of U /t. Doping rate 8 is different among 
(a)-(c): (a) in the overdoped regime, <5 ~ 0.20, (b) in the optimally doped and underdoped regime, 8 ~ 0.11, and (c) nearly half filling 
<5 ~ 0.03. In (d), a magnification of (c) in the vertical axis is shown to emphasize subtle variation of . 
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Fig. 29. (Color online) Momentum distribution function of d-wscve state for various values of t'/t for U/t = 12 and L = 12; the 
behavior basically does not alter if U/t ( ^ 10) or L varies. The doping rate is varied among the panels. The symbols for t'/t instructed 
in (a), (e) and (f) are common to all panels. Owing to the boundary condition, k y for each k-point shifts by n/(2L). The arrows with + 
[— ] in (c) and (d) indicate the direction of variation for t'/t > [t'/t < 0] when \t' /t\ increases; this t'/t dependence does not alter for 
(b), (e) and (f). The (7r,0) point is marked with a vertical guide line, and the (7r,7r) point with a small arrow on the abscissa. 
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around the antinodal point seems relevant. 44 ' 123 ' The 
maximum of |Vn(k)| is often used as an index of a quasi 
Fermi surface, and is reduced to a pure Fermi surface for 
|Vn(k)| — > oo. Here, we estimate |Vn(k)| near k = (tt, 0) 
(abbreviated as |Vn(X)|) from the finite-size VMC data 
as follows: dn(k)/dk x and dn(k)/dk y are obtained from 
finite differences of n(k) between X' and X given in §4.3, 
and those between X and X" = 7r(l, -j-), respectively. 

In Fig. 31, we plot n(k) for various 5's (L = 10-16) for 
three values of t'/t at U/t = 12. In contrast to the t'/t 
dependence [Fig. 29], an appreciable variation exists at 
any value of k, but the magnitude of variation in the anti- 
nodal area is still relatively large. We estimate |Vn(X)| 
for each S from Figs. 31(a)-(c), and plot it as a function 
of S with a triangle in Figs. 31(d)-(f), respectively, where 
5 dependence of P%° is displayed with circles together. 
For t' — 0, the behavior of |Vn(X)| closely agrees with 
that of P/f except near half filling, where, however, the 
disagreement obviously stems from the system-size de- 
pendence. 130 ' Thus, |Vn(X)| is almost proportional to 
P°° for t' = 0. For t'/t = 0.2, the behavior of |Vn(X)| 
is also nearly proportional to that of P//°\ |Vn(X)| ex- 
hibits similar system-size dependence near half filling. 
For t'/t = —0.25, the overall tendency of |Vn(X)| coin- 
cides with that of PS°. 



0.03 




U/t 



Fig. 32. (Color online) Comparison between SC correlation 
function P?° and Vn(k)| at k = X' ~ (tt, 0) as function of U/t for 
t' = and 8 = 0.102. The scales of the two quantities are adjusted 
so as to roughly equalize the maximum magnitudes. 



Finally, we touch on the U/t dependence of |Vn(X)|. 
If the simple ci-wave BCS function without correlation 
factors is used, we have |Vn(7r, 0)| = by differentiating 



i(k) = M 2 = 



£k - C 



( £k -cr+A 2 



(28) 



so that finite |Vrt(X)| is considered as a correlation effect. 
Figure 32 compares |Vn(X)| and P d °° of H d Q for t' = as 
a function of U/t in a slightly underdoped electron den- 
sity. The two quantities are both almost zero for U/t ^ 4, 
and broadly agree for any U/t, suggesting that the elec- 
tron distribution near (7r, 0) is important for c?-wave SC 
irrespective of the correlation strength. 

In this section, we have demonstrated that the 
strength of SC, Pj", behaves similarly to |Vn(X)| near 



regarding the stability of the d-wave state in previous 
sections, and is consistent with recent experiments, 32 ' 93 ' 
which argue that the antinodal electrons contribute to 
the SC gap. On the other hand, this result is contradic- 
tory to the concept of the dichotomy of electronic roles 
in the wave-number space, in which electrons near the 
nodal direction contribute to SC, and electrons near the 
antinodal point exclusively to the pseudogaps. 23 ' 24 ' Sim- 
ilar results were recently obtained by extended methods 
of DMFT. 131 ' 132 ' 

6. Summary 

With the high-T c cuprates in mind, we studied U/t, 6 
and t'/t dependence of a correlated c? 2; 2„j ) 2-wave singlet 
state and an antiferromagnetic (AF) state \&q F in- 
dependently as doped Mott insulators, applying them to 
a two-dimensional Hubbard (t-t'-U) model. The expec- 
tation values are calculated using a variational Monte 
Carlo method without additional approximations. Main 
results are recapitulated below. 

(1) As U/t increases, undergoes a sharp crossover 
of SC properties from a conventional BCS type to a 
kinetic-energy-driven type at U = U co somewhat larger 
than the band width (8t). As S decreases, U co continu- 
ously connected to the Mott transition point U c at half 
filling. For U < U co , superconductivity (SC), which is 
fragile, is enhanced by the pair scattering of the vector 
Q = (tt, tt) near the antinodes, according to the ordinary 
d-wave BCS theory. For U/t ^ 5, steady SC correspond- 
ing to the high-T c cuprates is not found. The SC corre- 
lation function P£° abruptly increases for J> 6. The 
unconventional SC for U > U co is robust for a wide range 
of t'/t and 6, which coincides with the range where the 
doublon-holon binding correlation — the essence of Mott 
physics — is effective, and the picture of the t-J model is 
valid. This SC cannot be understood by the instability 
of the Fermi surface in the underlying normal state. 

(2) By comparing the 5 dependence of d-wave SC cor- 
relation function P^° with the dome shape of T c and 
the condensation energy experimentally observed in the 
cuprates, we showed that the effective value of U for 
the cuprates must be larger than at least the Mott crit- 
ical value U c (~ 7t). Furthermore, we found from the 
analysis of kinetic energy that there are two kinds of 
holons for U > U co ; a holon created during a doublon for- 
mation does not contribute to conduction, but a holon 
introduced by doping becomes a carrier. Consequently, 
the number of doped holes (electrons) becomes equal to 
the carrier number. This is a natural picture in the t-J 
model, and consistent with various experimental results 
of cuprates. 116,117 ' On the other hand, for U < U co , all 
the holons (electrons) work in the same way as carriers, 
which make an ordinary large Fermi surface. This Fermi- 
liquid feature is contradictory to that of the cuprates. 

(3) In view of the two gap problem, the present result 
for U > U c is interpreted as follows: A gap for the sin- 
glet pair formation Ad, which possibly corresponds to a 
pseudogap, is a decreasing function of 5. Meanwhile, the 
strength of SC represented by P^° (related to the SC gap) 
has a dome-shaped 8 dependence, and is written by the 
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Fig. 31. (Color online) (a)-(c) Momentum distribution function of for various values of S at U/t = 12. The value of t'/t is varied 
among the three panels. The data for L = 10-16 are simultaneously displayed. Detailed instruction for 5 is omitted, but circles denote 
5 = 0, upward triangles <5 = 0.02-0.04, squares <5 = 0.06-0.11, downward triangles 8 = 0.12-0.14, upward pentagons 8 = 0.15-0.165, 
diamonds 8 = 0.165-0.21, stars 8 = 0.22-0.25, hexagons 8 = 0.26-0.29, and downward pentagons 8 = 0.29-0.31. To adjust the scale in 
the section (0, 0)-(7r, 0)-(7t, tt), the locus of k points between (it, tt)-(0, 0) slightly shifts depending on L. (d)-(f ) Comparison between SC 
correlation function PJ" and |Vra(k)| near k = (7r,0) estimated from the data in (a)-(c) for the three t' /Vs. We adjust the scales of the 
two quantities so as to roughly equalize the maximum magnitudes. The small digits for 8 ~ in (d) and (e) indicate the system size L 
used for the corresponding data point. 



bility such as the carrier density and the quasiparticle 
renormalization factor Z, which are increasing functions 
of 5 for small S. This suppression of charge fluctuation is 
imposed by the Mott physics. This feature of the SC gap 
is consistent with various theories for the t-J model. 

(4) A proper negative (positive) t'/t term stabilizes 
the SC (AF) state and destabilizes the AF (SC) state. A 
case of a very small 5 and a quite large U ft is exceptional. 
Meanwhile, the AF state is intrinsically unstable for hole- 
doped cases (t'/t ^ 0), and phase separates to the AF 
state with the local electron density n = 1 and the SC 
state with n < 0.85. 

(5) For the stability of AF state, the nesting of Q is 
primarily important on the whole Fermi surface. On the 
other hand, crucial for SC is the pair scattering of Q 
in the antinodal area, where the underlying band is flat. 
To check this aspect, we showed that the SC correlation 
P%° is closely connected to the slope of momentum dis- 
tribution function |Vn(k)| near the antinodal point. This 
contradicts the often argued dichotomy of electronic roles 
in the k space, i.e. the electrons in the nodal (antinodal) 
area exclusively contribute to SC (pseudogap). 

Finally, we briefly mention some future problems, (i) It 



is important to include long-range hopping terms in the 
band dispersion and the band renormalization effects ow- 
ing to the electron correlation, (ii) We should check the 
coexistence of SC and AF orders and their mutual ex- 
clusion, (hi) In connection with the pseudogap problem, 
normal states with some symmetry breakings should be 
studied, (iv) Making the most use of recent progress in 
VMC techniques, 133 -' we have to explore more exquisite 
trial wave functions. 114, 134 - ) (vi) We should study the role 
of doublons and holons in the SC phase more in detail, 
for instance, to check a recent proposal. 135 ' 
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Appendix A: Hubbard model for cuprates 

To consider the Cu02 planes, a proper starting model 
is the d-p model composed of 0-2p and C\i-3d x 2 - y 2 or- 
bitals. The d-p model is reducible to more simple models 
which still caDture the essence of CuOt Dlanes. Since a 
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Table B-l. Absolute values of total energy of the d- wave state, 
\E/t\, are compared among the three D-H correlation factors for 
8 > 0, t'/t = and L = 10. The last digits may have some errors. 
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Fig. B-l. (a) Local configurations around a doublon are clas- 
sified according to the number of holons in the nearest-neighbor 
sites, (b) A similar classification around a holon. 



doped electron enters a C\i-2>d x 2_ y 2 orbital, the model is 
reduced to a Hubbard model on a square lattice by elim- 
inating the degree of freedom for the 0-2p orbitals. In 
contrast, a doped hole enters an 0-2p orbitals; a Zhang- 
Rice singlet is formed, through which the model for low 
dopant density is reduced to the t-J model on a square 
lattice with small J/t. 2 ^ This t-J model is connected to 
the Hubbard model with large U/t (= At/ J) through the 
strong-coupling expansion 136 ) within the approximation 
of neglecting the pair-hopping (three-site) terms. Thus, 
the Hubbard model on a square lattice is derived as an 
effective model for both in electron-doped (ED) and hole- 
doped (HD) cuprates. Note, however, that the effective 
value of U ft in the Hubbard model may differ between 
ED and HD cases owing to the distinct derivation paths, 
even if the actual values of U in the Cu-3d x 2 _ y 2 orbitals 
are identical between the two. 

Another vital element of the model, eq. (1), is that 
one can map a more-than- half-filled case (n > 1) to a 
less-than-half-filled case (n < 1) by taking advantage of 
the antisymmetric level structure with respect to n = 1. 
The band dispersion, eq. (2), for k = (ir — k y , tt — k x ), 
which is the symmetric point of k with respect to the 
Fermi surface at half filling for t' — (AF Brillouin zone 
boundary), is written as, 

e(k) = 2i(cos k x + cos k y ) — At' cos k x cos k y , (AT) 

where tildes are the marks for k. To satisfy e(k) = — e(k), 
we need the relations t = t and r = —t', which are rep- 
resented by a canonical (electron-hole) transformation, 

cl->e«**< cur, (A-2) 

with Q = (tt,tt). Under this transformation, the model 
eq. (1) is invariant except for the addition of a constant 
U(N S — N). Consequently, calculations for an ED system 
with electron density n = 1 + 6 (S > 0) and t' Jt (< 0) 
can be replaced by those for the less-than-half-filled case 
with 1 — 8 and \t'/t\, if we consider t'/t is negative for 
both ED and HD cuprates. 

Appendix B: Details of doublon- holon factors 

In this Appendix, we compare different forms of the 
D-H correlation factor Vq [eq.(5)] to show that the dif- 
ference in results among Q® [eq. (7)] and other forms of 
Qj is insignificant. Here, Vq between nearest-neighbor 
sites is mainly considered, because Vq between further 
sites only make a minor difference. 20 '* 
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Fig. B-2. (Color online) The total energy of *q with Vf n (solid 
symbols) is compared with that with 'P? (open symbols) as a func- 
tion of U/t for four <5's. At half filling, Q^* H is reduced to Q| (solid 
stars). 17 ' Tendency for other values oit'/t is similar. 



[eq. (6)] (Vq) have been used, 80,81,137 - 1 on account of 
the particle-hole or doublon-holon symmetry. For S > 0, 
however, this symmetry is broken, and the number of 
holons always exceeds that of doublons. Consequently, 
if V^ is used, configurations S2, ^3 and 64 illustrated in 
Fig. B-l seem to appear frequently as 8 increases. In fact, 
a recent study 123 ' using a detailed parameterization for 
Qj has revealed that such configurations rarely appear 
for large U/t, because they increase the interaction en- 
ergy Ud. Anyway, we consider another asymmetric form 
of Qj, which is a natural extension of Q^: 

Q° H = »ddj JJ(1 - h j+T ) + (i h hj JJ(1 - d j+T ), (B-l) 

r r 

where binding parameters for holons hi, and doublons 
fid are optimized independently. At half filling, is 
reduced to Qj (fih, — Md) to retrieve the symmetry. As 
S increases, /ih is expected to decrease much faster than 
fid because £2 comes to appear less frequently than 82 ; 
therefore Q® n approaches Q° (fi^ rapidly decreases) . As 
shown in Table B-l, 138 - 1 the variational energy of Q® H is 
improved on that of Qj, especially for large U/t. In the 
following, we compare asymmetric forms and Q® . 

In Fig. B-2, the total energy of the d-wave state is 
compared between Q® and Qf H for t'/t = —0.25. In 
Table B-2, we list the energy improvement by Qj* H on 
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Table B-2. Improvement of total energy by ^q H on \Pq in per- 
centage, (E° — _B DH )/|i? D | X 100, for various parameter values and 
L = 10, which satisfy the closed-shell condition. 
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Fig. B-3. Weight assignment by Qj and by to local con- 
figurations with two doublons and two holons typical for U ~ U co 
at half filling. All the sites in the background are singly occupied. 
In each case, only relative positions of doublons to holons are illus- 
trated. In (a), there are two 8\ and two ei; in (b), 82, 80 and two 
£1; and in (c), £2, £0 an d two <5i- 
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(Color online) Comparison of optimized D-H correla- 

(/id and /ih) for nearest- 
neighbor sites. 



Qf (AE = E B - £ DH ) for some values of 5 and t'/t as a 
function of U/t. In every case, the energy is improved by 
Qf H . Comparing AE/t for t'/t = -0.25 and 0, we find 
that AE only slightly depends on t'/t. For each 5, AE/t 
has a maximum near U co /t = 10, especially at half filling, 
Qp H (=<3f) becomes highly advantageous for U > U c . 
For a fixed U/t, AE/t decreases as 8 increases. 

We briefly consider the reason why is advanta- 

geous at half filling and for U ~ U co . Assume that only 
two doublons and two holons exist in a local electron con- 
figurations at half filling; typical ones frequently appear- 
ing for U ~ U co are shown in Fig. B-3. The Gutzwillcr 
factor assigns a common weight g 2 to each of these con- 
figurations. As for the D-H factors, "Pq H and Vq give the 
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Fig. B-5. (Color online) Comparison of d-wave pair correlation 
function between two D-H correlation factors, Vq^ and "Pq. The 
behavior of PJ" is similar for other values of t'/t. 



(b), respectively; 123 ) however, for an unfavorable config- 
uration (c), Vq gives 1, but "Pq H reduces the weight as 
1 — /z. Thus, becomes better than at 6 = and 
for U ~ U co . When holes are doped, configurations like 
£2 (also £3 and £4) appearing in Fig. B-3(c) are rapidly 
suppressed; accordingly, the difference between Q™ and 
Q^j narrows. This aspect is clearly seen in the optimized 
D-H binding factors /i in and fid and fih in Qf K , as 
shown in Fig. B-4. As S increases, /id becomes slightly 
small, but /ih considerably and rapidly decreases, com- 
pared with /i in denoted by open symbols. 

Because the optimized values of the other variational 
parameters are similar between the two cases of Vq 11 
and Vq (not shown), physical quantities, e.g. n(k), 5(q) 
and N(q), thereof becomes similar (not shown). Among 
them, the d-wave pairing correlation function exhibit a 
relatively large difference as shown in Fig. B-5; the max- 
imum of is located at a slightly smaller U/t for Q 1)11 . 
However, the difference remains within quantitative. 

In conclusion, the D-H correlation factor Vq possesses 
properties sufficiently close to those of an improved factor 
P§ H at least for 6 > 0. 

Appendix C: Pairing correlation function 

We explain the details of how we determine long- 
distance values (Pj°) of the d^.^-wave SC correlation 
function Pd(r) defined by eq. (20) with eq. (21), to avoid 
misunderstanding. If Prf(r) remains finite for |r| — > 00, 
an off-diagonal long-range order exists. For finite sys- 
tems, however, one must be careful to estimate the value 
of Pd(r) for |r| = 00. 

In Fig. CT, the behavior of P<j( r ) f° r th° noninter- 
acting case {U/t — 0) is compared among different sys- 
tem sizes. Here, the path of r is chosen as shown in 
Fig. C-2(c). Because SC does not occur for U/t — 0, -Pd(r) 
must vanish for large |r|. However, for small systems like 
L = 10 and 12, -P<j(r) still has appreciable magnitude 
at distant points, e.g. r = (L/2,0) and r = (L/2, L/2). 
Furthermore, the location of large magnitude of Pz(r) de- 
pends on the values of t'/t and 6. If we require that Pd(r) 
should vanish for every distant point in this scale, we 
need to use a system with L <; 30. Such tendency remains 
for fairly large values of U/t, as seen in Fig. C-2. Fortu- 
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Fig. CM. (Color online) System-size dependence of d-wave SC 
correlation function Pd(r) for noninteracting systems (U/t = 0). 
The path of r is illustrated in Fig. C-2(c). The values of t'/t and 
<5 are (a) -0.25 and ~ 0.13, (b) and - 0.20, and (c) 0.2 and 
~ 0.13, respectively. 8 is chosen so as to satisfy the closed-shell 
condition under the same boundary condition as that in the VMC 
calculations. The data are obtained using the analytic formula. 



Fig. C-2. (Color online) U/t dependence of d-wave SC correla- 
tion function Pd(r) f° r L = 14 and <5 = 0.102. The path of r is 
depicted in (c); is not on this path. The value of t'/t is (a) 
—0.25 (hole doped), (b) 0, and (c) 0.2 (electron doped). In these 
systems, U max /t = 8.0, 8.1 and 10.5 for t'/t = -0.25, and 0.2, 
respectively. In every case, the closed-shell condition is satisfied. 
Note that the scale of the vertical axis is larger than Fig. CT. 
The data are obtained using the analytic formula for U/t = 0, and 
VMC otherwise. 



which is next to the farthermost point, is substantially 
zero (less than 10~ 4 ) at U = regardless of the values of 
t'/t, 6 and L, and remains smaller than those for r ^ 
for small U/t. Hence, we employ Pd(r oo) as P^° for small 
U/t (0 < U < t/ max ), where <7 max indicates the value of 
U of the largest Pd(r) for a fixed t'/t. Thereby, the arti- 
ficial enhancement of Pd(r) owing to the finite sizes can 
be eliminated. 

For large U/t (U > J7 max ), P%° can be determined 
more confidently. U/t dependence of Pd(r) for fixed L 
and S is shown in Fig. C-2. Pd(r) exhibits finite-size fluc- 
tuations for U Ss C/ max , but becomes almost constant for 
|r| > 3 139 > regardless of t'/t and L for U ^ C/ max . This 
behavior is the same as that found for the t-J model, 51 ) 
and reflects the short-range nature of strong correlation. 
Thus, we use the average value of -Pd(r) for |r| > 3 as 
P/T for U > U max . 

In summary, as a proper indicator of d-wave SC order, 
we safely define the long-distance value of P<j(r) for finite 



systems by, 



Pd(roo) 



(U < Umax) 



P2° = < 1 



M 



(C-l) 



M>3 



where r m = (i/2 — l,L/2), and M is the number of 
vectors Yj satisfying |r 3 -| > 3. 
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